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We use the canonical Hamiltonian formalism to generalize to spinning point particles the first 
law of mechanics established for binary systems of non-spinning point masses moving on circular 
orbits [Le Tiec, Blanchet, and Whiting, Phys. Rev. D 85, 064039 (2012)]. We find that the redshift 
observable of each particle is related in a very simple manner to the canonical Hamiltonian and, 
more generally, to a class of Fokker-type Hamiltonians. Our results are valid through linear order 
in the spin of each particle, but hold also for quadratic couplings between the spins of different 
particles. The knowledge of spin effects in the Hamiltonian allows us to compute spin-orbit terms in 
the redshift variable through 2.5PN order, for circular orbits and spins aligned or anti-aligned with 
the orbital angular momentum. To describe extended bodies such as black holes, we supplement 
the first law for spinning point-particle binaries with some "constitutive relations" that can be used 
for diagnosis of spin measurements in quasi-equilibrium initial data. 
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I. INTRODUCTION 

The need to develop faithful template waveforms for 
the search of gravitational waves from compact-object 
binary systems (compact binaries, for short) with cur- 
rent interferometric detectors [1, 2] has led, during the 
last years, to synergies and unexpected common ground 
among previously distinct research areas. The two main 
analytical frameworks used to study the rclativistic dy- 
namics of compact binaries are (i) the post-Newtonian 
(PN) approximation [3, 4], which describes the inspiral- 
ing motion beyond the Newtonian limit, in the small ve- 
locity and weak- field regime (» < c), and (ii) black hole 
perturbation theory [5], which describes both the weak- 
field and strong-field dynamics in the small mass-ratio 
limit (mi <C W2). The other main research area in the 
study of the two-body problem is numerical relativity 
[6], which aims at describing the strong- field regime by 
solving numerically the exact Einstein field equations. 
Furthermore, in order to build templates for the search 
of gravitational waves, the effective-one-body (EOB) ap- 
proach [7, 8] combines the information from those dif- 
ferent techniques in a flexible and effective way, and can 
provide accurate merger templates for advanced LIGO 
and Virgo. As said above, several studies [9-21] at the 
interface between those different frameworks have im- 
proved our knowledge of the two-body dynamics and 
gravitational- wave emission. 

Recently, Le Tiec, Blanchet, and Whiting [22] (hence- 
forth, Paper I) derived a "First Law of Mechanics" for 
binary systems of point particles moving on exact circu- 
lar orbits (compatible with an helical Killing symmetry). 
Using the first law, the authors found a very simple rela- 
tion between the PN binding energy of the binary and Dc- 
tweiler's redshift observable [11] — which can also be in- 
terpreted as the particle's Killing energy associated with 



the helical symmetry. This relation allowed the computa- 
tion of previously unknown high-order PN terms in the 
circular-orbit binding energy by using the redshift ob- 
servable computed numerically within the gravitational 
self-force (GSF) formalism [11-13, 23]. Moreover, the 
results in Paper I led to the following applications: 

1. Reference [19] computed the exact binding energy 
and angular momentum of two non-spinning com- 
pact objects (modeled as point particles) moving on 
circular orbits, at leading order beyond the test- 
particle approximation, and recovered the exact 
frequency shift of the Schwarzschild innermost sta- 
ble circular orbit induced by the conservative piece 
of the GSF [24]; 

2. Reference [20] built upon the works [14, 15, 19, 22] 
and derived the exact expressions of the EOB met- 
ric components through first order in the (symmet- 
ric) mass ratio. Quite interestingly, the results in 
Refs. [16, 18, 19] strongly suggest that the domain 
of validity of black hole perturbation theory calcu- 
lations may extend well beyond the extreme mass- 
ratio limit. 

Given the relevance of the first law of binary mechan- 
ics in enhancing our knowledge of the two-body dynamics 
in the comparable-mass case by using information from 
the perturbative GSF framework, we generalize it here 
to point particles carrying spins. We derive the first law 
of mechanics using the canonical Arnowitt-Deser-Misncr 
(ADM) formalism [25] applied to spinning point particles 
[26, 27]. In the PN context, both the ADM formalism 
and the harmonic-coordinate approach have been devel- 
oped up to high orders in the approximation, in order 
to compute the conservative part of the orbital dynam- 
ics of compact binaries. For non-spinning point masses, 
the Hamiltonian [28-30] and the harmonic-coordinates 



2 



equations of motion [31-35] are known through 3PN or- 
der. More recently, partial results at 4PN order have 
been reported [36, 37]. The logarithmic contributions to 
the conservative dynamics at 4PN and 5PN orders are 
also known [13, 22]. High-order spin-orbit effects have 
been computed in harmonic coordinates [38-40], in the 
Hamiltonian [26, 41, 42], and within the Effective Field 
Theory (EFT) approach [43]. See Ref. [44] for a review 
of the canonical ADM formalism for spinning point par- 
ticles, and its application in PN theory. In this paper we 
shall not consider spin-spin couplings or higher-order ef- 
fects in the spins. High-order spin-spin effects have been 
computed in the Hamiltonian [26, 45-48] and using EFT 
techniques [49-54]. 

The paper is organized as follows. In Sees. II and III 
we review the Lagrangian and Hamiltonian formalisms 
for systems of point particles carrying spins. In particu- 
lar, we derive a crucial relationship between the rcdshift 
observable and the variation of a general (Fokkcr-type) 
Hamiltonian, valid at linear order in the spins. In Sec. IV 
we derive the first law of mechanics for binary systems of 
spinning point particles using the canonical ADM formal- 
ism, and discuss some mathematical and physical con- 
sequences of this law. Then, in Sec. V, we employ the 
first law of mechanics and the ADM Hamiltonian to com- 
pute the spin-orbit contributions to the rcdshift observ- 
able through the (leading plus sub-leading) 2.5PN order. 
We find full agreement between our results and those re- 
cently obtained from a direct computation based on the 
near-zone PN metric [55]. In Sec. VI we discuss the first 
law in the particular case of binary black holes in corota- 
tion. Finally, Sec. VII summarizes our main results and 
discusses some prospects. Throughout this paper we use 
"geomctrized units" where G = c = 1. 



II. LAGRANGIAN AND HAMILTONIAN OF A 
SPINNING POINT PARTICLE 

In this Section, we review some necessary material for 
constructing a Lagrangian and then a Hamiltonian for a 
spinning point particle in curved spacetime. The formal- 
ism we shall use derives from the early works [56, 57]. It 
has recently been developed in the context of the EFT 
framework [49] . Alternatives and variants to this formal- 
ism can be found in Refs. [26, 44, 58]. The formalism 
yields for the equations of motion of spinning particles 
and the precession of the spins the classic results valid in 
general relativity [59-66]. 

Let us consider a single spinning point particle mov- 
ing in a given curved "background" metric g^v (x). The 
particle follows the worldline r M (s), with tangent four- 
velocity ~ dr^/ds, where s is a parameter along the 
representative worldline. In a first stage we do not require 
that the four- velocity be normalized, i.e., s need not be 
the proper time elapsed along the worldline. To describe 
some internal degrees of freedom to be associated with 
the particle's spin, we introduce a moving tetrad e a ll {s) 



along the trajectory. The tetrad is orthonormal, in the 
sense that g flu e ( /' e p = Va/3, and defines a "body- fixed" 
frame. 1 The rotation tensor f2 Ml< associated with that 
tetrad is defined by 

De 

-J^ = -We av , (2.1) 
ds 

where D/ds = u u V '„ is the covariant derivative with re- 
spect to the parameter s along the worldline. Equiva- 
lcntly, we have 

De v 

=e a »-^. (2.2) 
ds 

Because of the normalization of the tetrad the rotation 
tensor is antisymmetric: W = — fl vfi . 

A. Lagrangian formulation 

Following [49, 56], we look for an action for the spin- 
ning particle that is at once: (i) a covariant scalar, (ii) 
an internal Lorentz scalar, and (iii) reparametrization- 
invariant (i.e., its form must be independent of the pa- 
rameter used to follow the particle's worldline). We shall 
assume that the dynamical degrees of freedom are the 
particle's position and the tetrad ej 1 . We restrict 
ourselves to a Lagrangian depending only on the four- 
velocity m m , the rotation tensor f2 Ml/ , and the metric g^ v . 
In particular, this confines the formalism we are using 
to a polc-dipole model and to terms linear in the spins. 2 
Thus, the postulated "particle" action is of the type 

Vrt = /depart K,^,.g^) . (2.3) 

We cover the Lagrangian L part with a hat in order to 
distinguish it from the usual definition of the Lagrangian, 
ipart = ipart ds / dt (with t the coordinate time), that we 
shall use later. 

As it is written, depending only on Lorentz scalars, 
ipart is automatically a Lorentz scalar. By performing an 
infinitesimal coordinate transformation, one easily sees 
that the requirement that the Lagrangian be a covariant 
scalar specifies its dependence on the metric as being (sec 
e.g. Rcf. [57]) 

2 _^art = pV + gpnvp ( 24 ) 

og^u p 



Here T] a g = diag(— 1, 1, 1, 1) denotes the Minkowski metric. The 
indices fi, i>, ■ ■ ■ are the usual spacetime covariant indices, while 
a, /9, ■ • ■ are the internal Lorentz indices. The inverse (or dual) 
tetrad e a M , defined by ej 1 e a v = 8„ , satisfies rj a p e° (1 c'' ll = g^v 
We also have the completeness relation e ce M e' 3 = 
2 Such a model is "universal" in the sense that it can be used 
for black holes as well as neutrons stars. Indeed, the internal 
structure of the spinning body appears only at the next C(5 2 ), 
e.g. through the rotationally induced quadrupole moment. 
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where we have defined the conjugate linear momentum 
p^ and the antisymmetric spin tensor as 



dL 



part 



q _ <9-£<part 



aw 



(2.5a) 
(2.5b) 



Note that the right-hand side (RHS) of Eq. (2.4) is nec- 
essarily symmetric by exchange of the indices (i and v. 
Finally, imposing the invariance of the action (2.3) by 
reparametrization of the worldline, we find that the La- 
grangian must be a homogeneous function of degree one 
in the velocities and fl^ . Applying Euler's theorem 
to the function L part (M A ', fi pl/ ) immediately gives 



(2.6) 



where the functions p A i(it p , il pcr ) and S' A1I/ (u p , 51 P<T ) must 
be reparametrization invariant. Note that, at this stage, 
their explicit expressions are not known. They will be 
specified only once a spin supplementary condition (SSC) 
is imposed, as discussed in Sec. II C below. 

We now investigate the unconstrained variations of the 



action (2.3) with respect to ej' 



and g 1 *" '. First, we 



vary it with respect to the tetrad ej 1 while keeping the 
position r M fixed. We must have a way to distinguish 
intrinsic variations of the tetrad from those which are 
induced by a change of the metric g^ u . This is done by 
decomposing the variation 6e a v according to 



Se„" = 



<\fl 



1 



58^ + -5g 



(2.7) 



in which we have introduced the antisymmetric tensor 
56^ u = e a ^Se a "\ and where the corresponding symmet- 
ric part is simply given by the variation of the metric, 
i.e. e^^Sea = iSg^. Then we can consider the inde- 
pendent variations i56" il/ and 8g pv . Varying with respect 
to 69^ , but holding the metric fixed, gives the equation 
of spin precession which is found to be 



DS, 



d.S 



_ n pc _ O p <? 



pp > 



(2- 



or, alternatively, using the fact that the RHS of Eq. (2.4) 
is symmetric, 



D5 U 



(2.9) 



We next vary with respect to the particle's position r M 
while holding the tetrad e Q p fixed. Operationally, this 
means that we have to parallel-transport the tetrad along 
the displacement vector, i.e., we have to impose 



<$r"V„e c 



0. 



(2.10) 



A simple way to derive the result is to use locally inertial 
coordinates, such that T^ p = along the particle's world- 
line r p (s). Then, Eq. (2.10) yields SeJ 1 = 8r v d v e a ^ = 



— 8r v Y p p e a p = 0, and the variation gives the well-known 
Mathisson-Papapetrou [61-63] equation of motion 



Dj3p - 1 v »ft S*p ct 

j 2 fiVpG U 



(2.11) 



With more work, Eq. (2.11) can also be derived using an 
arbitrary coordinate system. 

Finally, varying with respect to the metric (while keep- 
ing S8 pl/ = 0) gives the stress-energy tensor of the spin- 
ning particle. We must take into account the scalarity of 
the action, as imposed by Eq. (2.4), and we obtain the 
standard result [59-66] 



rp t = J dsp^u 



§( 4 \x-r) 



V p / ds 5 p(p u v 



(2.12) 



where S^(x-r) denotes the four-dimensional Dirac func- 
tion, such that J d 4 x5^ 4 '(x) = 1. It can easily be checked 
that V„Tp!^ t = as a consequence of the equation of mo- 
tion (2.11) and the equation of spin precession (2.9). 



B. Hamiltonian formulation 

We now want to define a Hamiltonian associated with 
the Lagrangian (2.6). Because of the reparametrization- 
invariance condition, performing a Legendre transforma- 
tion with respect to the variables it p and Qf tv yields a 
vanishing result. Different routes are possible. One may 
add in the action a mass-shell constraint with a Lagrange 
multiplier and vary the action, keeping the momentum 
p M as an independent variable. Through various changes 
of variables and gauge fixings, the action is transformed 
into an action which possesses Euler-Lagrange equations 
of the form of Hamilton's equations, such that the canon- 
ical (ADM) Hamiltonian can be read off from these equa- 
tions (see Ref. [44] for a review). 

A different strategy, that we shall follow here, is closer 

to the usual procedure of classical mechanics. It consists 

of using as internal dynamical variables the six rotational 

degrees of freedom of Lorentz matrices, and defining the 

Hamiltonian by means of a 3 + 1 split, with a preferred 

time coordinate, with the usual Legendre transformation 

[58]. The tetrad e a p is decomposed into the product 

of an internal Lorentz transformation A' 3 and a refer- 
ci 

ence orthonormal tetrad field e a p {x), which is associated 
with the background metric g^ v {x) and is evaluated at 
the particle's position, x p = r M (s). 3 The point is that 
using as dynamical variables the six internal angles </> a 



3 Such decomposition was already implicit in the variation of the 
tetrad we performed to derive the spin precession equation (2.9). 
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of Lorcntz matrices allows us to have usual-looking 
Euler-Lagrange equations. We thus write 4 



,r)=A^)e/(r). 



(2.13) 



Next, we perform a 3+1 decomposition of all the fields. 
In particular, we split the particle's position and coordi- 
nate velocity according to r M = (i, r) and = (l,v), 
where v = dr/di and boldface letters denote ordinary 
spatial vectors (often also denoted with spatial indices i, 
j, ■ ■ ■ = 1, 2, 3). We choose for the parameter s along 
the worldline the proper time r and define the so-called 
rcdshift variable z by [11, 22] 



At 



I 1/2 



(2.14) 



Then the ordinary Lagrangian i par t = zL pa _ rtl such that 
•S'part = /dtLp ar t, becomes a function of the variables 
t, r, v, cj) a , and the ordinary derivative <p a = d(j) a /dt, 
with the dependence in t and r entering only through 
the background metric or, rather, through the reference 
tetrad evaluated at the particle's location: 



^part 



^part 



[v,0 o ,0 o ;e a "(t,r)] . 



(2.15) 



We have ordinary Euler-Lagrange equations for the gen- 
eralized coordinates (r, <f> a ) of the spinning particle (spa- 
tial vectors ranging over three components i = 1,2,3 
while the rotation label runs over a = 1, • • ■ ,6): 



(2.16a) 
(2.16b) 



dp t 


dL pait 


dt 


Q r l 




9L part 


dt 


d6 a 



where the conjugate momenta are defined in the usual 
way by Pi = dL pa , Yt /dv l and P^ a = dL palt /dij> a . They 
explicitly read 



<9A„ 



(2.17a) 
(2.17b) 



where 5*0,^ = e Q M e^ denote the tetrad components 
of the spin tensor in the "body-fixed" frame. The Euler- 
Lagrange equations (2.16a) and (2.16b) are equivalent 
to the equations of motion (2.11) and spin precession 
(2.9), respectively [58]. In particular, the spin precession 
equation is recovered in this formalism as the natural fact 
that the spin components remain constant in the frame 
attached to the body: 



dS a ^ 
~dV 



= 0. 



(2.18) 



4 The indices a, b = 1, ■ ■ ■ , 6 label the internal angles <p a of Lorentz 
matrices satisfying rj^g A 7 a A' 5 ^ = T) a f). Since e/ is a tetrad we 

have g^u e a M e^" = r} a p and, of course, still g^ u e p 



Va/3- 



Finally, the Hamiltonian is simply given in this ap- 
proach by the Legendre transformation as 



H, 



part 



P t v l + P 6 



v part j 



(2.19) 



and Hamilton's equations yield Eqs. (2.11) and (2.9), or 
equivalently (2.18) [58]. The Hamiltonian is a function 
of conjugate variables and again depends on the position 
r M only through the background tetrad: 



part 



-ffpart [Pi ■ 



,P<h 



(2.20) 



An explicit computation of the Hamiltonian (2.19) with 
the help of Eqs. (2.6) and (2.17) yields [58] 



part 



-pt 



l ( 

2-V 



Ouu t V tt a 



(2.21) 



Both terms should be understood as being functions of 
the conjugate variables, like in Eq. (2.20). As in ordi- 
nary classical mechanics, this is obtained by inverting 
Eqs. (2.17) to obtain v % and <j) a as functions of r l , Pi, </> a 
and P(f, a . However, here we need an additional relation 
to express pt as a function of the conjugate variables; this 
is provided by Eq. (2.6) which, in the 3 + 1 split, gives 



Pt 



L 



part ~ V l pi 



2 L 



C. Relating the Hamiltonian to the redshift 

In this Section we shall relate the Lagrangian -L part and 
the Hamiltonian H par t to the redshift variable z. To do so 
we must introduce a realistic physical model for the spin 
of the point particle in the polc-dipolc approximation. 

Up to now we have considered unconstrained variations 
of the action (2.3), which describes the particle's inter- 
nal degrees of freedom by the 6 independent components 
of the tetrad ej 1 — a 4 x 4 matrix subject to the 10 
constraints e.J^eJ' = r\ a p — or equivalently by the 6 
internal angles </> a . To correctly account for the number 
of degrees of freedom associated with the spin, we must 
impose 3 SSC. In this Section, we adopt the Tulczyjew 
covariant conditions [59, 601 



S» v p v = . 



(2.22) 



It would be possible to specify the Lagrangian in (2.3) so 
that the constraints (2.22) are directly the consequence of 
the equations derived from that Lagrangian [56]. Alter- 
natively, one could also introduce Lagrange multipliers 
into the action to enforce these constraints [44]. Here, 
for the sake of simplicity, we shall rather impose the con- 
straints (2.22) directly in the space of solutions. Fur- 
thermore, by contracting Eq. (2.9) with p v and using the 
equation of motion (2.11), one obtains the relation link- 



5 



ing the four-momentum p^ to the four-velocity as 5 



Ppipu) 



■m 2 u fl = -u x R\ p(r S'">S llv , 



(2.23) 



where we have defined m = —g pv p^p v . The parameter 
m is the mass of the particle, and it can be checked us- 
ing Eqs. (2.22) and (2.23) that it is constant along the 
trajectory, that is dm/dr = 0. 

Henceforth we restrict our attention to spin-orbit (SO) 
interactions, which are linear in the spins. Neglecting 
quadratic spin-spin (SS) and higher-order interactions, 
the linear momentum is simply proportional to the four- 
velocity: p^ = mu^ + 0(S 2 ). To linear order, Eq. (2.9) 
reduces to the equation of parallel transport for the spin 
tensor, DS^/dr = 0(S 2 ), and the corresponding stress- 
energy tensor reads 



rppV 

part 



V H V V 5(3) ( x 



+ 0(S 2 ) , (2.24) 



where 5^(x— r) is the three-dimensional Dirac function, 
such that J d 3 xS^ 3 \x) = 1, and we use the parametriza- 
tion by the proper time r and the notation (2.14), with 
v p = zu p . Furthermore, thanks to (2.6) the Lagrangian 
at linear order in the spin reads 



L 



part 



i 



0{S 2 



(2.25) 



We now want to compute from Eq. (2.25) the partial 
derivative of the Lagrangian with respect to the mass 
m of the particle, while holding the dynamical degrees 
of freedom r, v, <fr a , and </> a fixed. The first term in 
(2.25) will contribute in an obvious way; however we must 
also control the partial derivative of the spin tensor S^ V1 
which is the conjugate of the rotation tensor f^", with 
respect to the mass. Let us show that, in fact, the second 
term in Eq. (2.25) will not contribute to the result at the 
required O(S). 

To show this we choose the time-like tetrad vector to 
agree with the four- momentum rcscalcd by the mass, i.e. 
e fl = p p jm. Then, using Eq. (2.1) and the equation of 
motion (2.11) we get Wp v = \u v W v S pa . Comparing 



with the SSC (2.22) and using p^ 



mu,. 



-0(S 2 



we in- 



fer that Sp V and f2^„ must be proportional, up to a small 
curvature coupling to the spin and some higher-order spin 
terms: oc Q^„ - ■^Rp V p t7 S pcr + 0(S 3 ). For instance, 
such a proportionality relation is verified in models for 
the relativistic spherical top [49, 56], and the constant 



5 We denote (pu) = p v u v . By further contracting Eq. (2.23) with 
M M we obtain an explicit expression for (pu), which can then be 
substituted back into (2.23). 



of proportionality is associated with the moment of in- 
ertia, say /, of the spherical top. The previous relation 
can be solved iteratively, yielding an expression for Sp V 
as a function of m, /, R 



pupa > 



and Q pl ' , which schemati- 
cally reads S~l[tl- ^RCl + (j^f RRQ + ■■■]. The 
functional dependence on m is explicit in this expression. 
Taking a partial derivative with respect to the mass, and 
replacing back W in terms of the spin, we deduce that 
dSpv/dm = O(S). Since the rotation tensor is also O(S), 
we finally get WdS^/dm = 0(S 2 ). Wc thus conclude 
that the second term in Eq. (2.25) will not contribute at 
linear order in the spin. Thus, we have proven within the 
covariant SSC (2.22) that 



dL 



part 



dm 



-z + 0{S 2 ). 



(2.26) 



Using the properties of the Legendrc transformation 
(2.19), we readily find that the partial derivative of the 
Hamiltonian with respect to the mass, holding the con- 
jugate variables r, Pi, cf> a , and P^ a fixed, is 



dH, 



part 



dm 



z + 0{S 2 ). 



(2.27) 



The latter result will be central to the derivation of the 
first law of mechanics for binary systems of spinning par- 
ticles in Sec. IV. For the moment we note that it has 
been proven only for a single spinning particle moving in 
a given curved "background" spacetime. 



III. LAGRANGIAN AND HAMILTONIAN OF 
INTERACTING SPINNING PARTICLES 

In this Section, we consider a self-gravitating matter 
system consisting of N spinning point particles, which we 
shall label by A, B = 1, • • • ,N. The Lagrangian for each 
particle is of the form (2.15) in a given metric, so the sys- 
tem of N particles (interacting only through gravitation) 
is described by the matter Lagrangian 



A 



l A ,m A ;e a »(r A )]. (3.1) 



A=l 



Here we have added for later convenience the depen- 
dence on the mass = (— p^p'X) 1 ^ 2 - Neglecting terms 
quadratic or higher order in the spins, this Lagrangian 
reduces to 



In 



JV 

£ 

A=l 



~-A 



-m A 



2 P" A 



0(S 2 A ) 



(3.2) 



In the above equation, the spins S^ v should be seen as 
functions of r Al v A , (j) A , and <j) A through their relations 
to the rotation tensors as discussed in the previous 
Section. Next, we add the usual Einstcin-Hilbert term 
for gravitation, Leh = J <i 3 x^~gR, and obtain the 
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total Lagrangian describing the "matter + gravitation" 
system as 

L = L m&t \\ A ,<j)\,$\,m A \eJ l {r A )\ + L EK [g^„] , (3.3) 

in which the spacetime metric reads g^ v = e^e^, and 
the brackets refer, with obvious notations, to a func- 
tional dependence. Varying the Lagrangian (3.3) with 
respect to the tetrad ej 1 we get the Einstein field equa- 
tions G^ v = 87rT^ t , where the matter tensor is a sum 
of terms of the type (2.24), one for each particle. 



The Fokker Lagrangian [67] is now obtained by insert- 
ing the formal solution (3.5) into (3.4); hence 



t harm r 

-t-'F — ^mat 



[v A , 4>% 4>% m A ; eJ*(r A )] + L^ m [fa] 



(3.7) 



This Lagrangian gives the correct equations of motion 
(and spin precession) for the matter variables. Its vari- 
ational derivative with respect to either the position r A 
or spin parameters <f) A [which we collectively denote as 
cr A = (r A ,(j) A )] reads 



A. Fokker Lagrangian 

Wc shall define the Fokker Lagrangian [67] by for- 
mally solving the Einstein field equations and eliminat- 
ing the gravitational field degrees of freedom in Eq. (3.3). 



Writing the metric as g^ p 



the field equa- 



tions are solved perturbatively in powers of by us- 
ing the graviton propagator. As usual we discard in the 
Einstcin-Hilbcrt action a total derivative and consider 
instead the Landau-Lifshitz Lagrangian, which involves 
only first-order derivatives of the metric; symbolically 
Lll = J d 3 Xy/^gVr. Furthermore, we need to spec- 
ify a coordinate system to solve the field equations; this 
is done by adding some gauge-fixing terms to the action. 
Choosing the harmonic coordinate system, this yields the 
"harmonic-gauge-relaxed" Lagrangian -L^l™ which con- 
tains at quadratic order the graviton propagator V^pa-, 
say Llf' m = -\ J d 3 x [hV^h + 0(h 3 )} (sec e.g. App. A 
of Rcf. [68]). We thus consider 



Lharm j 
— i'mat 



[va,« 



> A ,(p A ,m A ]e 



- harm 
7 LL 



(3.4) 



Next, we perturbatively solve the harmonic-gauge Ein- 
stein field equations and we formally rc-cxpand the solu- 
tion in a PN approximation. Expanding the retardations 
of retarded integrals in a PN scheme yields accelerations 
and time derivatives of accelerations. It will also generate 
higher time derivatives of the spin variables. The solu- 
tion, valid at any point x M in the near-zone, is a function 
of all the source parameters, 

^( x ) = e a »[x-,r B ,v B ,!iB,(t> a B,<i> a B,$ a B,rriB] , (3.5) 

where we collectively denote the higher time derivatives 
by &b = (dvg/df , d 2 v#/di 2 , • • • ) and symbolically by <f> B 
for the spin variables. Then the near-zone metric solution 
is guv = tap, £ a v - Thus, by construction the solution (3.5) 
satisfies 



5<T A 



5er A 

6e» SL harm 



VB,<f> B ,<pB, m B;e a V{r B )] 



6(T A 5€a 



(3.8) 



where the last term vanishes thanks to Eq. (3. 6). 6 The 
equations of motion and precession are precisely those we 
expect for spinning particles within the background e a M 
generated by all the particles themselves, i.e. 



6L* 



8 a A 



■[v B ,^,^,m B ;e/(r B )] =0. 



(3.9) 



To obtain this result, the Fokker Lagrangian (3.7) must 
involve not only the matter part evaluated in the back- 
ground generated by the particles, but also, crucially, the 
Einstcin-Hilbcrt (or Landau-Lifshitz) part for the gravi- 
tational field. 

We can now use the properties of the Fokker La- 
grangian to compute its partial derivative with respect 
to one of the masses, say m A , holding all dynamical vari- 
ables (and the other masses) fixed. As before we obtain 



5jL harm 

drriA 



dm a 
dej 1 SL haim 
drriA <5e Q M 



vs,0B,0B,Wij;e a M (rij)] 



(3.10) 



where again the last term vanishes thanks to Eq. (3.6), 
such that it remains only 



aL harm 

drriA 



8L 



part 



dm,A 



v i 



(t> A ,j> A ,m A ;e a i*(r A )], (3.11) 



where we used Eq. (3.1) to express the result in terms of 
the single-particle Lagrangian for particle A. From (3.2) 
and the proof given above Eq. (2.26), we have, at linear 
order in the spin S A , 



aL harm 

dniA 



-za + 0{S\). 



(3.12) 



8L 



harm 



5L n 



Se c 



■[v A ,(j)%^X,mA;eJ l (rA)] 
[fa] = ■ 



6 Notice that in the last term of Eq. (3.8) the functional derivative 
(3.6) fita^/dcA is to be taken in a generalized sense, because of the 

presence of accelerations ag and (fig. 
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We emphasize that this result is valid in the background 
field generated by the system of N particles; the redshift 
variable therein reads in full glory 

/ ■ •• \V 2 

Z A = (-g^[r A ;rB,VB,a.B,(t>%^%,<p a B,m B ]v'Xv A ) 

(3.13) 

In particular, it depends on the masses ms of the other 
particles in the system. 

B. Fokker Hamiltonian 

That a relationship analogous to Eq. (3.12) holds also 
for the Hamiltonian is not completely straightforward, 
because of the presence of higher-order derivatives (no- 
tably accelerations) in the harmonic-coordinate Fokker 
Lagrangian (3.7). The higher-order derivatives come 
from the replacement of the metric by the perturbative 
solution (3.5), and from explicit time derivatives in the 
self-interaction terms (such as dh dh) of the graviton ac- 
tion. However, as shown in Ref. [69] (see also [33, 70]), 
one can remove from a PN expanded Lagrangian the 
higher-order derivatives of through a suitable redefi- 
nition of the position variables, r A — > r^ ew , and the ad- 
dition of a total time derivative to the Lagrangian. This 
process can be implemented to all orders in PN theory 
[33, 69, 70]. The perturbative method of Ref. [69] is very 
general, and there should be no obstacle in applying it 
to the spin variables, i.e. in removing from a Lagrangian 
the higher-order derivatives of 4> a A through a redefinition 
4> A —> (4 l A ) ncw and the addition of a total time derivative. 
Thus we can eliminate all the higher-order time deriva- 
tives in the harmonic-coordinate Fokker Lagrangian (3.7) 
by introducing the new (non-harmonic) Lagrangian 

r T harm £ j harm J s~i 

t rharm , oij F r . 01j F X ia , aU fn i A \ 

L F = L F + -^—8v A + -^54> A + — , (3.14) 

where we denote 5r A = y a cw — v A and 8<j> A = {(f> A ) nevr — 
(j> A , and G is some function of time that does not affect 
the dynamics. The Lagrangian (3.14) now depends on 
r A ew and v A ew , but no longer on accelerations & A ew ; and it 
depends on the spins only through (4> A ) nev ' and (<j> A ) new - 
Also, the non-harmonic metric (or tetrad) solution of the 
field equations in the new variables takes the same form 
as in Eq. (3.5), but without the contributions from higher 
time derivatives. This is because the redefinition of the 
positions r A — > r A cw can be seen as being induced by a 
coordinate transformation of the "bulk" near-zone met- 
ric, and that metric when evaluated at the location of 
the particles is necessarily free of higher derivatives (see 
Sec. 3 in Ref. [70]). Thus, we have the same result as in 
Eq. (3.12). We shall from now on remove the superscript 



As usual, non-linear contributions in accelerations can be re- 
moved by the addition of so-called "multiple-zero" terms [71, 72], 



"new" on all variables, but we shall keep in mind that 
they correspond to non-harmonic coordinates. 

Finally, we have an ordinary Lagrangian L-p depend- 
ing on the dynamical variables (r A ,v A ,(f) A ,(f> A ); like in 
Eq. (2.19) we apply a standard Legendre transformation 
to obtain the corresponding Hamiltonian, say H-p. Defin- 
ing = dL-p/dv l A and = dLp/d^) A , we have 

N 

H ? = E { p ? v ' 1 a + pfA) - l f > ( 3 - 15 ) 

A=l 

where v A and 4> A are now viewed as functions of the 
conjugate variables rg, Pf, <f) a B , and P?, obtained by 
inverting the definitions of the momenta P t A and P^ . 
And like in Eq. (2.27), a consequence of the Legendre 
transformation (3.15) is that the partial derivative of the 
Hamiltonian with respect to the mass m A , holding the 
conjugate variables rg, Pf , </>^, and P? (and the masses 
tub 7^ Tn A ) fixed, is 

^= z A + 0(S 2 A ). (3.16) 
am A 

The result (3.16) is the main ingredient of the first law of 
mechanics for binary systems of spinning point particles, 
which is the topic of Sec. IV. 



C. Constructing the canonical Hamiltonian 

Working within the ADM Hamiltonian formalism, the 
authors of Refs. [26, 27] showed that, starting from the 
unconstrained Hamiltonian depending on the spin vari- 
ables S^, it is possible to build the constrained canonical 
Hamiltonian for gravitationally interacting particles by 
imposing suitable constraints and gauge conditions, and 
making a change of variables for the positions, momenta, 
and spins of the particles. On the other hand, starting 
from the Hamiltonian (2.19) which depends on the spin 
variables <fi a and P<p a , the construction of the constrained 
canonical Hamiltonian for a single spinning particle in a 
given curved spacetime was explicitly done in Ref. [58], 
thereby extending to curved spacetime the classic work 
by Hanson and Regge [56] for a relativistic top in flat 
spacetime. Being guaranteed that a constrained, canoni- 
cal ADM Hamiltonian exists for gravitationally interact- 
ing spinning point particles [26, 27, 44], henceforth we 
assume that the construction of Ref. [58] can be extended 
to interacting particles. Let us then briefly review how 
the constrained Hamiltonian for a single spinning particle 
was built in Ref. [58]. 

We start from the unconstrained Hamiltonian for a sin- 
gle spinning particle H pai t(r,Pi,(f> a ,P^ a ), together with 
the Poisson bracket operation {• , •} (see Sec. Ill B in 
Ref. [58] for the explicit expressions of the Poisson brack- 
ets between the relevant variables). We consider a set of 
constraints £fc = with k = ,2M (with M < 9) 
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such that the matrix Cki = is not singular. Fol- 

lowing Ref. [56], these constraints can be imposed by (i) 
replacing the original Poisson brackets with Dirac brack- 
ets, 8 and (ii) inserting the constraints directly in the 
original, unconstrained Hamiltonian. More specifically, 
Ref. [58] first imposed the Newton- Wigner SSC [56] , suit- 
ably generalized to curved spacetime, to eliminate 3 out 
of the 6 variables P<f, a (a = 1, ■ ■ • ,6), and an additional 
constraint to eliminate 3 out of the 6 variables 4> a , such 
that the constraint hypersurface contains the same num- 
ber of configuration coordinates and conjugate momenta. 
Then, the constrained Hamiltonian Hc(r,p,S) and the 
phase-space algebra of the constrained system were com- 
puted, showing that they are canonical at linear order in 
the particle's spin. Thus, the variables r, p, and S — with 
S i = y ijk S jk where S jk = efefS^ 1,2,3)- 
obey the standard commutation relations {r l ,pj} = 5j 
and {S^S^} = e^ k S k , all other brackets vanishing. As 
usual, this algebra can be used to compute the dynami- 
cal evolution of any (time-independent) function / of the 
canonical variables according to / = d//di = {/, H}. Fi- 
nally, the constrained canonical Hamiltonian constructed 
in Ref. [58] coincides with the ADM Hamiltonian for a 
spinning particle moving in a given background. 

As said, we shall assume that the construction above 
can be extended to gravitationally interacting spinning 
point particles. The resulting constrained Hamiltonian 
Hc{ya, Pa, Sa) with A = 1, • • • ,N with usual canonical 
variables will coincide with the ADM canonical Hamil- 
tonian. Furthermore, being constructed from the class 
of Fokker Hamiltonians (3.15), it necessarily satisfies 
Eq. (3.16). Thereafter, we restrict ourselves to a binary 
system of spinning point masses, i.e. A = 1,2, and we 
denote the canonical ADM Hamiltonian by H. 



IV. FIRST LAW OF MECHANICS FOR BINARY 
POINT PARTICLES WITH SPINS 

A. Derivation of the first law 

In Section IIIB we proved the remarkable relation- 
ship (3.16), valid when neglecting non- linear terms in the 



spins of the type O(Sf) and 0(S$ 



i.e. 



dH 

dniA 



ZA ■ 



(4.1) 



Note that this relation is not only valid for linear contri- 
butions 0(5*1) an d 0(6*2), but also for non-linear spin 
contributions of the type ©(SiS^). Here the rcdshift 



For two phase-space functions A and B, the Dirac brackets are 
given by {A,B} D = {A, B} + {A, £fc}{-B, £i}[C _1 ] w , and repre- 
sent the projection of the original symplectic structure onto the 
phase-space surface defined by the constraints. 



quantity za is defined in terms of the coordinate four- 
velocity v^ A = (1, \a) of particle A, in ADM-type coordi- 
nates, by 



ZA 



1/2 



(4.2) 

with being the background metric generated by the 
two particles, and evaluated at the location of particle A. 
Note that the particle velocity va in the RHS of Eq. (4.2) 
should be viewed as a function of the conjugate variables 
ta, Pa, Sa, obtained by inverting the definition of the 
momentum pa- We emphasize that Eqs. (4.1)-(4.2) are 
valid for generic orbits and spin configurations. For cir- 
cular orbits and spins aligned or anti-aligned with respect 
to the orbital angular momentum, za would coincide with 
the (gauge-invariant) redshift observable introduced by 
Detweiler [11]. This is the redshift of a photon emitted 
from the particle and observed at a large distance, along 
the symmetry axis z perpendicular to the orbital plane. 

In order to derive the first law, we reduce the ADM 
Hamiltonian H{va, Pa, Sa, wa), given in a generic frame 
of reference, to the center-of-mass Hamiltonian, say 
H (r, p, Sa, tua)- This is achieved by first imposing that 
the ADM linear momentum Padm = Pi + P2 vanishes, 
then substituting the individual momenta pa for the rel- 
ative linear momentum p = pi = — P2, and denoting the 
coordinate separation by r = — r 2 . 

We shall also limit our study to spins aligned or anti- 
aligned with the orbital angular momentum L, and to 
circular orbits. Defining L = iz, where z is the unit vec- 
tor orthogonal to the orbital plane, we have Sa = Sa z, 
with I Sa I < ni>A- Using polar coordinates (r,tp) in the 
orbital plane, we can express the Hamiltonian as a func- 
tion of the separation r = |r|, of the radial momentum p r , 
of the azimuthal momentum p v = L, and of the masses 
toa and amplitudes Sa of the spins of the two point par- 
ticles. An unconstrained variation of the Hamiltonian 
H(r,p r ,p v ,SA,mA) thus gives 



STT dH . dH . 
oH = — — dr + —— dp r 

or op r 

dH . 

orriA + 



dH 
dH 



5L 



\dmA 



dS A 



SS A 



(4.3) 



where to ease the notation we have omitted to indicate 
that the partial derivatives with respect to r, p r , L, 
Sa, and toa are computed keeping all the other vari- 
ables fixed. If the variation compares one solution of the 
Hamiltonian dynamics to a neighboring solution, then 
Hamilton's equations of motion must be satisfied. For 
circular orbits, this yields dH/dr = —p r = 0, as well as 
dH/dp r = f = 0. The constant circular-orbit frequency 
is given by dH/dL = <p = f2, and the Hamiltonian is nu- 
merically equal to the ADM mass: H = M. Therefore, 
"on shell" we have 



^ ( dH . dH 

22{d^ 6m - + dSA- 



SS A 



(4.4) 



9 



At linear order in the spins, the Hamiltonian can be 
written as the sum of an orbital part Hq , which does not 
depend on the spins Sa, and a spin part H$ = J2a ^aSa- 
The partial derivatives of the Hamiltonian with respect 
to the spins, 



dH 
dSA~' 



(4.5) 



are the so-called precession frequencies of the spins. In- 
deed, using the algebra satisfied by the canonical vari- 
ables, the spins Sa can easily be seen to satisfy — in the 
most general, precessing case — the Newtonian-looking 
(but exact) precession equations Sa = ft a X Sa- Thus, 
the usual Euclidean norms of the canonical spins are con- 
served, Sa ■ Sa = const, and Sa are also referred to as 
the constant-in-magnitude spins, or constant spins. 9 

Finally, combining Eqs. (4.1), (4.4), and (4.5), the first 
law of mechanics takes on the simple form 1 



5M — flSL = (z A Sm A + ft a SS a ) 



(4.6) 



This differential relation gives the changes in the ADM 
mass M and orbital angular momentum L of the binary 
system under small changes in the individual masses tua 
and spins Sa of the two point particles. We emphasize 
that Eq. (4.6) is valid for any spin magnitude |£a| < m A . 

It is convenient to replace the orbital angular momen- 
tum L in favor of the total (ADM-like) angular momen- 
tum J. For aligned or anti-aligned spins, we simply have 



J = L 



5>. 



(4.7) 



In terms of J, the first law (4.6) reads 

5M -Q6J = Y^UASm A + {^A-^)SS A ■ (4.8) 

A 

We recall that the expressions (4.6) and (4.8) of the first 
law are, in principle, valid only at linear order in the am- 
plitudes Sa of the constant (or canonical) spin variables 
Sa- It is true that the definitions (4.5) of the precession 
frequencies do not require that the Hamiltonian H be 
linear in the spins; if contributions quadratic (or higher) 
in the spins were included in the Hamiltonian, then the 
precession frequencies £Ia would become functions of the 
spins Sa- Nevertheless, our proof in Sec. Ill of the result 
(3.16) is only valid at linear order in the spins. 



9 Note that the constant spins are not uniquely defined, because a 
local Euclidean rotation leaves the magnitude ■ Sa of the spin 
Sa unchanged. However, for spins aligned or anti-aligned with 
the orbital angular momentum, this remaining gauge freedom is 
irrelevant, as it does not affect the (algebraic) magnitude Sa- 
10 A similar result was previously established in [73], see Eq. (24) 
there, although without the crucial mass variation terms za Stua 
in the right-hand side. 



B. Consequences of the first law 

We now explore some consequences of the first law for 
spinning point particles. In particular, we establish alge- 
braic expressions which can be regarded as first integrals 
associated with the variational relations (4.6) and (4.8). 

For two spinning point particles with masses tua and 
spins Sa, on a circular orbit with azimuthal frequency 
ft, the ADM mass M, orbital angular momentum L, red- 
shifts za and precession frequencies ft a arc all functions 
of the five independent variables (ft, mi, mi, Si, Si)- 
Therefore, the first law (4.6) is equivalent to the following 
set of partial differential equations: 



— -n— 
~dft~ dfi 



dM 

diriA 
dM 



ft 



ft 



dL 

dniA 
dL 



dS A dS A 



0. 



ZA , 



(4.9a) 
(4.9b) 
(4.9c) 



Note that the relations (4.9a) and (4.9b) arc unchanged 
with respect to the results in the non-spinning case [see 
Eqs. (2.40) and (2.41) of Paper I], where they have been 
checked up to 3PN order plus the logarithmic contribu- 
tions at 4PN and 5PN order. However, in the spinning 
case we have the additional equations (4.9c) relating the 
precession frequencies ft a to the partial derivatives of M 
and L with respect to the spins Sa- 

By combining Eqs. (4.9a)-(4.9b) with numerical cal- 
culations of the GSF effect on the redshift z\ of a non- 
spinning point particle on a circular orbit around a 
Schwarzschild black hole [11, 23, 74], the exact expres- 
sions of the total mass M and orbital angular momentum 
L could be determined, at leading order beyond the test- 
particle approximation [19]. This result was used in [20] 
(see also Ref. [21]) to improve the knowledge of the EOB 
model [7, 8, 75] for non-spinning binaries. We leave to fu- 
ture work the generalization of these findings to spinning 
point masses. 

By commutation of partial derivatives, Eqs. (4.9b) and 
(4.9c) can be combined to give the interesting relation 



3za Ofls 



dS E 



dniA 



(4.10) 



This equation relates the variation in the redshift za of 
particle A under a small change in the spin Sb of particle 
B to the variation of the precession frequency fts of par- 
ticle B under a small change in the mass tua of particle 
A. Therefore, for A 7^ B, Eq. (4.10) reflects some equi- 
librium state of the spinning point particles under their 
mutual gravitational attraction. For A = B, we obtain 
a non-trivial relation between the redshift (coinciding, in 
adapted coordinates, with the helical Killing energy) and 
the precession frequency of each spinning point particle. 

We now derive the first integral associated with the 
variational first law (4.6) by following the same steps as 
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in the proof given in the non-spinning case in Paper I. 
First, we introduce the convenient combination 



M = M - QL . 



(4.11) 



Making the change of variables (Q, mi,ni2,Si, S2) — > 
(f2, m, v, xii X2)i where m = m\ + m% is the total mass, 
v = mi 7712/ 'm the symmetric mass ratio, and \A = 
Sa/iti\ the dimensionless spins, Eqs. (4.9b)-(4.9c) can 
be combined to give 



^(waza + 2CIaSa) = ( m A 



dM 

1 dm a 



dM 

dm 



(4.12) 



Next, we notice that the ratio M/m is dimensionless and 
symmetric by the exchange mi f> 7772 of the particles; it 
must thus be a function of mfl, v, xii an d X2] see e.g. 
Eqs. (5.7) below. This last observation implies the rela- 
tionship md(M/m)/dm = Q,d(M/m)/dCl, which when 
combined with (4.9a) and (4.12) yields the first integral 
relation 



M 



2QL = ^2{m A z A + 2n A S A ) 



(4.13) 



Alternatively, Euler's theorem for homogeneous functions 
provides a more straightforward proof of the result (4.13). 
Since the Einstein field equations do not contain any priv- 
ileged mass scale, the ADM mass must be a homogeneous 
function of degree one in L 1 / 2 , 777,4 and S 1 / 2 . Hence, Eu- 
ler's theorem implies 



M(L,m A ,S A )=L 1 / 2 



+ 



A 



dM 

dL^ 
dM 

m A - 



' dm a 



S 



1/2 



dM 



dS 



1/2 



(4.14) 



which when combined with the first law (4.6) immedi- 
ately yields the first integral relation (4.13). This result 
generalizes to spinning point particles the first integral 
relation derived in the non-spinning case [see Eq. (1.2) of 
Paper I], as it involves appropriate additional spin terms. 

Furthermore, we can use the previous results to obtain 
a relationship between the spin contributions to the ADM 
mass M and the sum J2a 171 aza of the redshifted masses. 
Indeed, Eqs. (4.10) and (4.9c) successively imply 
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dn f 



dS, 



A 



m A z A 



^ mA dm a 



dS, 



777A- 



dM 



dmA 
(4.15) 

By combining Eqs. (4.12) and (4.13), we can easily re- 
place the sum J^a mAdM/dmA in the RHS by the ex- 
pression -M + 2M - 2 J^a SAdM/dSA- At linear order 
in the spins, 9 2 A^/(95a<95's) = 0, and we find 



d 
dS^ 



(71/ + 777lZl + 777lZ2) = . 



(4.16) 



This simple relation immediately gives the SO contribu- 
tions to the ADM mass M once those in the redshift 
obscrvables Z\ and zi are known, and vice versa. 

Until now, all results were given in terms of the or- 
bital angular momentum L. We can rewrite them using 
the total angular momentum (4.7) instead. The partial 
differential equations (4.9) then become 



dM 
~dQ 



dM 

dmA 
dM 



0^ 

dn 

dJ 

dmA 



0, 



ZA , 



' ) 7 



dS A dS A 



(4.17a) 
(4.17b) 
(4.17c) 



Equation (4.17a) is the "thermodynamical" relation com- 
monly used in PN theory for quasi-circular orbits (see, 
e.g., Refs. [29, 76]), or in the construction of sequences 
of quasi-equilibrium initial data for binary black holes 
and binary neutrons stars [77-80]. In terms of the total 
angular momentum, Eq. (4.13) becomes 



n 



m - 2nj = y m A z A - 2(n - n A )s A 



(4.18) 



which is, as expected, the first integral associated with 
the variational first law (4.8). The existence of such a 
simple, linear, algebraic relation between the local quan- 
tities za and Qa on one hand, and the global quantities 
M and J on the other hand, is noticeable. 



V. SPIN-ORBIT EFFECTS IN THE PARTICLE'S 
REDSHIFT OBSERVABLE 

We employ the ADM Hamiltonian to derive the spin- 
orbit terms at ncxt-to-leading 2.5PN order in the parti- 
cle's redshift observable (or helical Killing energy) using 
the relation (4.1). We keep neglecting all terms 0{S\) 
or higher. We do not include non-linear spin interactions 
of the type ©(SiS^) either, even though they must sat- 
isfy our first law of mechanics. The center-of-mass ADM 
Hamiltonian can then be written as 



H = H(r,p,SA,m A ) = Ho+H i 



so 



NLSO 



(5.1) 



where Ho is the orbital (or non-spinning) Hamiltonian, 
which is known through 3PN order, while the leading- 
order 1.5PN spin-orbit term reads [81] 



(5.2) 



with L = r x p the orbital angular momentum, n = r/r 
the unit vector pointing from 7772 to mi, and 



S eff = ( 1 



37772 
47771 



Si 



1 



37771 
47772 



(5.3) 
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The next-to-leading 2.5PN SO terms, computed first in 
the harmonic-coordinates approach [38, 39], and then in 
the ADM Hamiltonian [41], read 11 



2.5PN order. To check those relations we use the spin- 
orbit contributions to the binding energy E = M — m 
and orbital angular momentum: 



-^nlso — 



(L-S X ) 



-6mi — 13m2 



3 J_ 

4 m 2 
3 1 
\~m\ 

1 « — > 2. 



1 



2 mi TO2 
7 1 



4 mi m 2 



(n p 

5 mi 
8mf 



mi 

2 



(5.4) 



As was done in Sec. IV, we now restrict ourselves to 
circular orbits (p r = n • p = 0, fl = dH/dL, dH/dr = 0) 
and spins aligned or anti-aligned with the orbital angular 
momentum, and use Eq. (4.1) to compute the spin-orbit 
coupling through 2.5PN order in the redshift observable 
za- Remember that the partial derivative with respect to 
m,A in Eq. (4.1) is to be taken at fixed r, L, Sa, and ms 
tua- A straightforward calculation gives the following 
leading 1.5PN SO contributions to the redshift associated 
with particle l: 12 



v so 
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(5.5) 



and similarly for the next-to-leading 2.5PN SO terms 
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(5.7a) 
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(5.7c) 



7 



7 , 847 217 , 

-A -I v A v \i 

2 144 144 72 



-I-Ia 
2 2 



847 
144' 



2!7 A 7 2 

Av v 2 

144 72 



X2 



,5/2 



(5.7d) 



We need also the orbital (i.e. spin-independent) part of 
the precession frequencies at next-to-leading order. We 
compute it from the definition (4.5) together with the SO 
Hamiltonian (5.2)-(5.4), and find 
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We have introduced the usual frequency-related PN pa- 
rameter x = (mf2) 2 / 3 and we recall the notation xa = 
S A /m 2 A . The SO and NLSO terms found in Eqs. (5.5) 
and (5.6) agree with the results of Ref. [55] , obtained from 
a direct calculation based on the near-zone PN metric; see 
Eq. (3.13). This provides an important, non-trivial test 
of the validity of the relation (3.16) derived in Sec. III. 

For consistency we verified that Eqs. (4.9), (4.10), and 
(4.13) hold in the spin-orbit sector. We also verified that 
Eq. (4.16) is satisfied by the spin-orbit terms through 
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Note that when computing the partial derivatives to 
prove Eqs. (4.9), (4.10), and (4.16), we need to hold fixed 
the spin variables Sa and not the reduced spins xa- 

We conclude that the redshift obscrvables za satisfy 
the first law of mechanics in the non-spinning sector, up 
to 3PN order included, as well as for the leading-order 
4PN and next-to-leading order 5PN logarithmic terms 
(Sec. II D of Paper I), and, in the spin-orbit sector, up 
to the next-to-leading 2.5PN order. 



VI. BINARY BLACK HOLES IN COROTATION 



11 Note that spin-orbit terms in the ADM Hamiltonian are also 
known at 3.5PN order [42, 48], but we do not include them here. 

12 Hereafter we assume, without any loss of generality, m\ ^ m,2- 
We denote the reduced mass difference by A = (m2 — mi)/m = 
VI — 4ia The redshift observable 22 of particle 2 can immedi- 
ately be deduced from z\ by setting A — > — A and xi 4 — ► X2- 



In this Section, we consider the particular case of coro- 
tating point particles on a circular orbit. By modeling 
each spinning point particle by an isolated rotating black 
hole, and comparing with the first law of binary black 
hole mechanics [82], we define and compute the proper 
rotation frequencies of the corotating point particles. 
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A. Spinning point particles as Kerr black holes 

In Sec. IV we derived the first law of mechanics for 
spinning point particles [see Eq. (4.8)]. We now derive 
an alternative version of this result that can heuristically 
be applied to binary black holes carrying spins. Clearly, 
in order to derive a first law for extended bodies such 
as rotating black holes, we must supplement the pole- 
dipole model we have used so far with some "constitutive 
relations" myt(/M; Sa,- ■ • ) specifying the energy content 
of the bodies, i.e. the relations of their masses tua to 
the spins Sa and to some "irreducible" masses fiA- In 
principle, these constitutive relations should also involve 
other parameters such as some external (tidal) multipolc 
moments due to the environment; in the two-body case, 
they would thus depend on the orbital separation and 
the mass of the other body. 

By analogy with the first law of mechanics for a single 
black hole [83] , we define for each spinning point particle 
the analogue of an irreducible mass \ia = rn l A via the 
variational relation 5mA = caS^a + waSSa, in which 
the "response coefficients" ca and uja are associated with 
the internal structure: 



ca 



dmA 



dfj,A 



UJA 



dmA 



dS A 



(6.1) 



In particular, lo a can be interpreted as the proper rota- 
tion frequency of body A. Assuming the validity of the 
Christodoulou mass formula m A = fx A + S\/(A^l a ) for 
Kerr black holes [84, 85], i.e. neglecting the influence of 
the other body, we have 



ca 



U1 A 



HA 
TUA 

S A 



1 



°A 

4/4 



4mAH 2 A 



(6.2a) 
(6.2b) 



For an isolated black hole, the coefficient ca is related to 
the constant surface gravity ka by ca = ^Ha k a- Note 
that, at linear order in the spins, we have ca = 1+0(S A ). 
While the first law (4.8) does not account for spin ef- 
fects 0(S A ) or higher, we shall not make the substitution 
ca — > 1 in the equations and discussion below, in order 
to ease the comparison with the binary black hole case. 

Identifying each spinning point mass with an "ex- 
tended body" (or black hole) characterized by coefficients 
ca and uja, we use the definitions (6.1) and find that the 
first law (4.8) can be rewritten as 



5M-n8J 



A 



\caza S^a + {za + Oa - f2) SSa 



(6.3) 

This variational relationship is reminiscent of the first 
law of binary black hole mechanics of Friedman, Uryii, 
and Shibata [82]: 



5 M — 8 J = y j 4/ia k a Sua 
a 



(6.4) 



which holds for two actual black holes with irreducible 
masses fiA and constant surface gravities ka (the surface 
areas of the black holes being A a 



167T/4) 



B. First law for corotating binaries 

Both first laws (6.3) and (6.4) have been derived for 
circular orbits and spins aligned with the orbital angu- 
lar momentum. However, whereas Eq. (6.3) is valid for 
arbitrary spin magnitudes, Eq. (6.4) can only describe 
corotating black holes. This key difference is intimately 
related to the fact that black holes are finite-sized objects, 
whereas point particles have (by definition) no spatial ex- 
tension. The assumption of a circular orbit implies the 
existence of a helical symmetry. In the binary black hole 
case, the helical Killing field must be tangent to the null 
geodesic generators of the horizons, which entirely con- 
strains the rotational state of each black hole. If corota- 
tion were not realized, the resulting non- vanishing shear 
would lead to the growth of the horizon's surface areas, 
in contradiction with the hypothesis of helical symmetry 
[82]. In the binary point-particle case, however, no such 
restriction occurs because of the point-like nature of such 
idealized objects. 

By analogy with Eq. (6.4), we sec from Eq. (6.3) that 
the two point particles will be "corotating" if and only if 



ZAOJA + Qj 



n. 



(6.5) 



This equation determines the values of the proper fre- 
quencies uja to be used to compute the corotating point 
particle spins Sa through the relation Sa = ^tiaHa^a 
(and the Christodoulou mass formula). Physically, the 
condition (6.5) means that the redshifted proper rotation 
frequency of each particle, zauja, must be equal to the 
circular-orbit frequency f2, as seen in a frame rotating at 
the angular rate fi^ with respect to an inertial frame of 
reference. Alternatively, we can associate to the proper 
rotation frequency uja a proper rotation angle <f>A by the 
definition uja = d^x/drA, such that Eq. (6.5) becomes 
d(f>A = (fi — &a) dt — dtp — d&A- This relationship gives 
the change d<f>A in the proper rotation angle in terms of 
the changes dip and d$^ in the orbital phase and preces- 
sion angle during a coordinate time interval dt. 

For such corotating point particles, the first law (6.3) 
simplifies considerably to 



SM — 0,5 J = caza Sua 



(6.6) 



It is almost identical to that derived for non-spinning 
binaries in Paper I, which is Eq. (6.6) with the substi- 
tutions 5 ha SrriA and ca — > 1. Since the irreducible 
mass fj-A of a rotating black hole is the spin- independent 
part of its total mass m,A, this observation suggests that 
corotating binaries are very similar to non-spinning bi- 
naries, at least from the perspective of the first law of 
mechanics. This helps to explain why the first law of bi- 
nary black hole mechanics (6.4) has been observed — in 
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the context of sequences of quasi-cquilibrium initial data 
— to be remarkably well satisfied by non-spinning binary 
black holes [80], even though it should, in principle, hold 
only for corotating black hole binaries [82]. 

We also note that the variational first law (6.6) admits 
the first integral 

M - 20J = ^2 VACaza , (6.7) 

A 

which again is very similar to the first integral relation de- 
rived in the non-spinning case, i.e. Eq. (6.7) with (ia —> 
rriA and Ca — > 1- Equation (6.7) can also be derived from 
the first integral (4.18) valid for arbitrary spins, in which 
we replace the masses tua by the irreducible masses \xa 
and spins Sa using ttla = ^aca + ^aSa (the analogue of 
Smarr's formula [86] for spinning point particles treated 
as Kerr black holes), and impose the condition (6.5) for 
corotation. 

Now, comparing the first law (6.4) for corotating black 
holes with the first law (6.6) for corotating point parti- 
cles, we notice the formal analogy 

c a za < — ► 4/j, a ka ■ (6.8) 

This point-particle/black- hole analogy was pointed out in 
Paper I (with ca — > 1), in the case of non-spinning point 
masses. Up to the factors of ca = 1 + 0(S\). the first 
law for corotating point particles expressed in terms of 
the irreducible masses ha is identical to the first law for 
non-spinning point particles expressed in terms of the to- 
tal masses itia = ha + 0(S\). As derived in Paper I, the 
analogy (6.8), with ca —> 1, was not entirely physically 
motivated, because it required the identification of the 
irreducible masses of the corotating black holes with the 
total masses of the non-spinning point particles. Having 
been established here in the case of corotating point par- 
ticles, and requiring only the identification of irreducible 
masses, the relation (6.8) is much more compelling. 

In the limit of large separation, the black holes or point 
particles can be viewed, in first approximation, as iso- 
lated. In that limit, we know that A^a^a ca for each 
black hole, and za — > 1 for each point particle, such that 
the analogy (6.8) is consistent. Going beyond the large 
separation limit, Eq. (6.8) suggests that the deviations 
of the redshifts of the point particles from one provide a 
measure of the interaction between the black holes. In 
particular, since the redshifts za = dr^/dt must be less 
than one, we expect 



in binary systems of corotating black holes on circular 
orbits. In other words, the tidal interaction between 
the holes should decrease the surface gravities with re- 
spect to their values in isolation. This prediction could 
be checked by computing numerically ka in sequences 
of quasi-circular initial data relying on the existence of 
a global helical Killing vector [77, 78]. One could also 



compare quantitatively PN results for za [11-13, 22] with 
numerical results for A/j,aka, as functions of the circular- 
orbit frequency fl. 

C. Proper rotation frequencies through 2PN order 

The condition (6.5) for corotation can be solved for the 
proper rotation frequencies uja of the particles: 

uja = u\ (n - A ) , (6.10) 

where u\ = 1/za has been computed up to very high PN 
orders, for non-spinning binaries [11-13, 22]. Here we 
only need the lPN-accurate results. For circular orbits, 
and in the center-of-mass frame, the spin- independent 
contributions to the redshift associated with particle 1 
read 

uS = 1+(| + |a-0* (6.11) 
/27 27 A 5 5 A v 2 \ 9 3 , 

However, for corotating binaries we should also include 
spin contributions. The lowest order spin-related effect 
comes from the leading-order 1.5PN SO terms, as given 
by Eq. (5.5). Since \a oc f2 at leading order, this term 
will yield a 3PN contribution, which can be neglected 
here. Combining Eqs. (5.8) and (6.11), the condition 
(6.10) for corotation readily fixes at 2PN order 

w A = O jl - vx + v + x 2 + 0(a; 3 ) j . (6.12) 

Noticeably, even though u\ ^ u\ and f^i ^ fl2, we find 
uii = u>2 up to 2PN order included. Future work should 
investigate whether this symmetry property still holds at 
the next 3PN order. Note that in the Newtonian limit 
x — > or the test-particle limit v — > we simply have 
uja = in agreement with physical intuition. 

The observation that u>\ = ui2 up to high order suggests 
that the proper rotation frequency uja might physically 
correspond to the rotation rate of the tidal field of the 
companion B ^ A, as measured in the local asymptotic 
rest frame (LARF) of body A [87]. By matching the near- 
zone PN metric of a binary system of point particles to 
the metric of two tidally distorted Schwarzschild black 
holes, Alvi [88] (see also Ref. [89]) computed the 1PN- 
accurate expression of the rotation rate of the tidal field 
of body B in the LARF of body A. Interestingly, his 
result agrees with the restriction at IPN order of our 
2PN-accurate expression (6.12). 

Furthermore, Caudill et al. [80] previously made that 
same hypothesis, in the context of quasi-cquilibrium ini- 
tial data for equal-mass black hole binaries. They found 
that using the lPN-accurate result uja = £1(1 — vx) in- 
stead of the leading-order result uja = £1 considerably 
improved the agreement between two different quasi-local 
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measures of the individual spins of the black holes; see 
Eqs. (54)-(56) and Figs. 4 and 5 of Ref. [80]. It would 
be interesting to revisit their work using our improved. 
2PN-accurate expression (6.12) for the proper rotation 
frequencies. 

In Ref. [76] , the binding energy E = M — m and total 
angular momentum J of a binary system of corotating 
point particles were computed by replacing the masses 
niA and spins Sa by the irreducible masses fj,A and 
proper rotation frequencies Ua, using the PN-expanded 
form of the relations between (m,A, Sa) and (ha, wa) ob- 
tained from the Christodoulou formula, and assuming the 
leading-order result oja = 0- Using the corrected for- 
mula (6.12) instead, the results of Ref. [76] are modified. 
Redoing the calculation for the additional spin-related 
contributions E cor and J cor to the binding energy and 
total angular momentum in the corotating case, we find 

6Tj)a£+ W (-10 + 2577)4, 



M (2 



(6.13a) 



fi J cor = At (4 - 12??) xl + fj,ri (-16 + 40r?) a£ , (6.13b) 

where the total mass /j, = fj,± + fj,2 , the symmetric mass ra- 
tio <r\ = //i/i2/ 'm 2 \ and the dimensionlcss invariant PN pa- 
rameter = (/i£l) 2 / 3 are now expressed in terms of the 
irreducible masses [ia, rather than the masses m^- The 
2PN and 3PN spin- related contributions (6.13) must be 
added to the known 3PN-accurate expressions of E and J 
for spinless binaries. As expected, the 1PN correction in 
Eq. (6.12) modifies the 3PN terms in (6.13) with respect 
to the results of Ref. [76] . It can be checked that the ad- 
ditional contributions (6.13) for corotating binaries now 
satisfy the relation 



dE c 



O 



d,J c 



on 



(6.14) 



Since the spin-independent contributions to E and J are 
known to satisfy that same relation (see Paper I), we 
conclude that the thermodynamical law is verified, as it 
should be according to the first law (6.6), for corotat- 
ing point particles. It would be interesting to revisit the 
PN/NR comparison of Ref. [76] using the corrected for- 
mulas (6.13). 



VII. SUMMARY AND PROSPECTS 



Similarly to previous work [19, 20], a calculation of the 
gravitational self-force effect on the rcdshift observable of 
a particle on a circular, equatorial orbit around a Kerr 
black hole [55, 90] could be used to compute the exact 
spin-orbit contributions to the binding energy and orbital 
angular momentum, for aligned or anti-aligncd spins, be- 
yond the test-particle approximation. This would allow 
computing the frequency shift of the Kerr innermost sta- 
ble circular orbit under the effect of the conservative self- 
force (at linear order in the spin), an important strong- 
field benchmark. 

New comparisons to numerical relativity simulations 
of spinning black hole binaries could then explore fur- 
ther the promise of using perturbation theory to model 
comparable-mass compact binaries [16, 19]. This infor- 
mation could also be used to improve the EOB model 
for spinning binaries [91-96]. Furthermore, once second- 
order gravitational self-force calculations [97-100] be- 
come mature enough, these same ideas could be applied 
to compute the fully relativistic second-order contribu- 
tions in the (symmetric) mass ratio to the binding energy 
and angular momentum, for circular orbits. All of these 
results would be valuable to improve template waveforms 
for inspiraling and coalescing binaries of compact objects. 

Moreover, our simple relation between the redshift ob- 
servable and the (ADM, or more generally Fokker-type) 
Hamiltonian allowed us to compute spin-orbit effects in 
the redshift through 2.5PN order, for circular orbits and 
spins aligned or anti-aligned. These results agree with 
a recent direct calculation based on the PN metric [55]. 
Future work can extend the knowledge of the rcdshift 
function to higher PN order using the recently computed 
spin-orbit effects at 3.5PN order [40, 42]. It would also 
be interesting and useful to extend the first law to sec- 
ond order in the spins, in order to account for spin-spin 
contributions and quadrupolar deformations. 

Finally, we specified our first law for binary systems of 
spinning point particles to the corotating case, allowing 
a comparison with the binary black hole case. Extending 
previous results [88, 89], we computed the proper rota- 
tion frequencies of the particles through 2PN order. This 
finding could be used to improve the quasi-local measure 
of the individual spins of black holes in the context of 
quasi-equilibrium initial data [80] . 



In this paper we generalized to spinning point particles 
the first law of binary mechanics established in Paper I 
for non-spinning point masses. We derived a simple rela- 
tion between the Hamiltonian and the redshift observable 
using a general formalism for the Fokkcr Lagrangian and 
Hamiltonian of a system of spinning particles. We then 
derived the first law within the canonical ADM Hamilto- 
nian formalism, for binaries on circular orbits and with 
spins aligned or anti-aligned with the orbital angular mo- 
mentum. We also obtained several useful relations link- 
ing the main quantities describing the two-body dynam- 
ics for circular orbits, mainly the ADM mass, angular 
momentum, and precession frequencies. 
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